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Abstract. We compute the anomalous dimensions of the flavour non-singlet twist-2 Wilson 
operators in the RI'/SMOM scheme at two loops in an arbitrary linear covariant gauge. In 
addition we provide the full Green's function for these operators inserted in a quark 2-point 
function at the symmetric subtraction point. The three loop anomalous dimensions in the 
Landau gauge are also derived. 
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1 Introduction. 



In a quantum field theory the behaviour of the Green's functions or n-point functions derived 
from the Lagrangian carry all the information about the dynamics of the quantum particles. For 
the vast majority of quantum field theories, however, it is impossible to extract their behaviour 
for all ranges of momenta and parameters, such as the particle masses and coupling constants. 
Instead one invariably examines them in various regions of interest, such as at high or low en- 
ergy, using a variety of techniques. For instance, at high energy in Quantum Chromodynamics 
(QCD) one uses perturbation theory since the coupling constant is small as a consequence of 
asymptotic freedom, [U [2]. By contrast at low energy perturbation theory breaks down and 
non-perturbative methods have been developed and refined to give credible information. The 
central tool at this energy is lattice gauge theory involving an intense amount of numerical com- 
putations on high performance computers. This approach has in general been hugely successful 
in determining bound state masses, for example, and exploring the structure of nucleons. For 
instance, matrix elements of the underlying operators used in deep inelastic scattering, known 
as Wilson operators, play a key role in this, [3]. The behaviour of such matrix elements at low 
energy is useful in extracting theoretical information for parton structure functions. In outlining 
the general aspects of Green's functions in understanding the dynamics of the strong nuclear 
force, we are overlooking the huge technical effort which is required to ensure accurate estimates 
are obtained. For instance, as one is dealing with a quantum field theory, the operators undergo 
renormalization, [3]. Equally when one produces estimates from a low energy computation on 
the lattice one has to be assured that the result is consistent with and extrapolates onto the 
high energy behaviour of the same object which can be computed within perturbation theory. 
Indeed there has been a large degree of activity on the lattice in this respect for quark currents 
and Wilson operators. For instance, seefilSElElElElIinillllliaiiaiiaiiailSlIITlIIHlI^ 
for representative analyses. 

For the continuum computations one usually calculates in the MS scheme which is a mass 
independent renormalization scheme. The advantage of this scheme is that it is one in which 
the largest order of perturbation theory can be determined. However, defining the same scheme 
for lattice calculations is not as easy since it invariably requires a numerical differentiation on 
the lattice. Taking such derivatives carries a financial penalty. Therefore, alternative schemes 
have been developed for lattice gauge theory which avoids the use of derivatives. Such a class 
of schemes is generally referred to as Regularization Invariant (RI), [2H I22j . However, there 
are two main types. The original class involves RI itself and a modified version known as RI', 
|2H I22j. These are similar in that QCD is renormalized for 3-point and higher Green's functions 
according to the MS prescription but for 2-point functions such as those determining propagator 
corrections, the renormalization condition is defined by ensuring that the contributions from 
radiative corrections at the subtraction point are absent. In this class of schemes we include the 
zero momentum insertion of an operator in, say, a quark 2-point function. The modified scheme, 
RI', differs from RI in the way the quark wave function is defined, [21^ I22j . Whilst originally 
defined in \21\ [22] specifically in the context of the lattice, this scheme has also been studied 
in the Landau gauge in the continuum at three and four loops, [23], and at three loops in an 
arbitrary linear covariant gauge, [21]. Subsequently, the Green's functions of a variety of low 
moment twist-2 flavour non-singlet operators central to deep inelastic scattering inserted in a 
quark 2-point function were evaluated to three loops in MS and RI', [251 126] • This high order 
of perturbation theory provided useful information on matching the lattice measurement of the 
same object at high energy. 

More recently a second class of regularization invariant schemes has been developed in [271 
[281 [29]. It is termed RI'/SMOM where the first part of the designation indicates the RI' scheme 
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definition of the quark wave function renormalization. The second refers to the method used 
to define the renormahzation constants of 3-point functions with an operator insertion at non- 
zero momentum. As there are momenta flowing into each point of the Green's function, the 
renormahzation is carried out at a symmetric point where the squared momenta of ah three 
incoming momenta take the same value which is the origin of the S. The MOM indicates the 
ethos that was present in the overall RI definition in that the condition for defining the operator 
renormalization constant is to ensure that the part corresponding to radiative corrections at 
this symmetric subtraction point are absent. This particular scheme was developed to avoid the 
strong sensitivity of the RI' scheme to infrared effects, [27], as well as to try and have a more 
rapidly convergent perturbation series for the conversion functions to the MS scheme. It was 
initially applied to the renormalization of the quark currents at one loop in |27] to determine the 
anomalous dimensions, amplitudes and conversion functions in the MS scheme. Subsequently, 
the two loop renormalization of the scalar and tensor currents was treated in [28\ [29] with the 
anomalous dimensions and conversion functions being computed. Though with the knowledge 
of the two loop conversion functions the three loop Landau gauge anomalous dimensions were 
deduced too. More recently, the full set of amplitudes for the scalar, vector and tensor currents 
have been calculated at two loops in [30]. 

Given that the measurement of nucleon matrix elements is of interest it is the purpose of 
this article to provide the anomalous dimensions and amplitudes for the second moment of 
the flavour non-singlet Wilson operator at two loops in RI'/SMOM. This will build on the 
analogous one loop computation of [31]. The treatment of this operator is complicated by the 
fact that it mixes with a total derivative operator. Ordinarily one is not concerned by this 
extra operator since if it were inserted at zero momentum it would give no contribution to the 
Green's function. However, the symmetric subtraction point of the RI'/SMOM scheme means 
that such total derivative operators play an active role and cannot be ignored. The full three 
loop MS mixing matrix of anomalous dimensions for this second moment operator was given 
in [32j. With the two loop conversion functions we will provide the three loop result in the 
RI'/SMOM scheme. One feature of the RI'/SMOM construction for the Wilson operators is the 
role the vector current has in the renormalization. Its relation to the Slavnov- Taylor identity was 
discussed in [27] at one loop, as well as two loop, [30], and in the context of the Wilson operators 
in [31j. The divergence of the vector current resides within the total derivative operator into 
which the Wilson operator mixes. Hence, the two loop renormalization of the vector current 
and the associated RI'/SMOM amplitudes of the Green's function where this was determined, 
|30j . are important as checks on the computation presented here. Indeed partly related to this 
is that one has to ensure the renormalization of the second moment is consistent with one of the 
Slavnov- Taylor identities of QCD. Further, the third moment of the Wilson operator was also 
considered at one loop in |31] . From that it transpires that the renormalization of the higher 
moment operators is dependent on that of the lower moment operators of the tower with the 
vector current at the foundation. Therefore, the treatment of the n = 2 moment is relevant for 
that of the higher moment operators. Throughout we will use flavours of massless quarks and 
therefore all our expressions are in the chiral limit. Including masses for quarks is not currently 
viable as the basic scalar master Feynman integrals have not been evaluated for the momentum 
configuration of the Green's function we consider here. 

The article is organized as follows. The background to the formalism we will use, notation 
and method to define and extract the scalar amplitudes of interest are discussed in section two. 
Section three is devoted to summarizing the main results of the RI'/SMOM computation of the 
second moment of the Wilson operator. The amplitudes are given numerically. Though the 
full exact expressions are recorded in the form of a set of tables. The results relating to the 
actual operator renormalization, such as the anomalous dimensions and conversion functions are 
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provided in section four whilst we summarize our conclusions in section five. 



2 Preliminaries. 

First, we recall the key properties of the operators we will consider. We use the notation of ^32j 
in this respect and denote 

W2 = Sipj^'D"^ , dW2 = Sdi^ii^Y^) (2-1) 

where Z)^ is the covariant derivative with all derivatives acting to the right and S indicates that 
the operator is traceless and symmetric in its Lorentz indices. Moreover, as we will be using 
dimensional regularization we regard the operators as being traceless in d-dimensions. So, for 
instance, we have 

SO^,^ = O'^J + O^^^ - ^v^.O^''' (2.2) 

where 

O^,' = i^^^D,^ . (2.3) 

The total derivative operator denoted by dW2 also obeys the same symmetrization. Throughout 
we concentrate on flavour non-singlet operators and omit the flavour indices on the quark fields 
tp. Ordinarily when one renormalizes W2 in perturbation theory the operator is inserted in 
a quark 2-point function, [3| 1331 [34], at zero momentum. Then the renormalization constant 
emerges under the assumption that the renormalization of this operator is multiplicative. This 
is not strictly speaking true. It is only true for that particular momentum routing through the 
quark 2-point Green's function with the operator insertion. If instead there was a momentum 
flowing out through the operator then there is a mixing with the total derivative operator. The 
Feynman rule for the latter operator vanishes when there is no momentum flow out through the 
operator which is why the momentum configuration of [31I331IM] is multiplicative. In other words 
for the full renormalization of W2 one has to determine the mixing matrix of renormalization 
constants where dW2 completes the set. As we are considering the renormalization of W2 in the 
RI'/SMOM scheme we have to take this into account as the momentum setup for the underlying 
Green's function involves two independent momenta flowing in through the external quark legs 
and out through the operator itself. This is illustrated in Figure 1. 



t p+q 




Figure 1: The Green's function yp{p)0''^^ {—p — q)tp{q) 

The renormalization of the set of level W2 operators has been carried out to three loops in the 
MS scheme in [32]. As we will require these later for the three loop RI'/SMOM scheme anomalous 
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dimensions later we briefly recall the relevant results and formalism. First, irrespective of which 
scheme one works in the mixing matrix of renormalization constants is 



Z. 



W2 







^12 
7VK2 



(2.4) 



Here we retain the notation of [3H [32] whereby the superscript indicates the operator level W2 
and the labels indicating the matrix elements respectively label the two elements of the operators 
in that set which are W2 and dW2- Whilst W2 is used to mean level and operator it will be clear 
from the context which is meant. The anomalous dimensions are formally defined by 



(2.5) 



where here and below W2 represents the level, from which one can deduce the relations between 
the individual matrix elements are 







7iT^(a,aX^ 

7i?(a,a)Z,t 
72^^(a,a)Z2^^ 



d_ 



d 



+ 7l2 (O, ajZ22 + ^^^12 



-Z. 



W2 
22 



with 



P{a)— + a-/a{a,a)- 



(2.6) 



(2.7) 



dfi ' ^ 3a da 
and 7fj(a, a) is the anomalous dimension of the gauge parameter. We retain our previous con- 
ventions, [311 [32], for consistency. In the MS scheme since the operators W2 and dW2 are gauge 
invariant then the anomalous dimensions are independent of the gauge parameter. However, for 
mass dependent renormalization schemes, of which RI'/SMOM is an example, the anomalous 
dimensions are a dependent. Therefore, we recall that for the MS scheme 



7it(a) 



1 



MS 



-Cpa + — [376CaCf 



112C7? 



+ 



243 



(5184C(3) + 20920) ClCp - (15552C(3) + 8528) CaC], 



(10368C(3) + 6256) CACpTrNf + (10368C(3) - 560) 



(10368C(3) - 6824) ClTpNf - SmCpT^N^ 



7l^^(a) 



MS 3 



-Cpa + — 



27 



56C7| - 188CaCf + GiCFTpNf 



+ 



243 



(7776C(3) + 4264) CaC^ - (2592C(3) + 10460) C^Cf 
+ (5184C(3) + 3128) CaCfTpN} - (5184C(3) - 280) 0% 
- (5184C(3) - 3412) C^TpNf + USCpT^N^] + 0{d 



MS 



(2.8) 



where C{z) is the Riemann zeta function and for this scheme we have omitted the part of the 
argument involving a to emphasise the absence of gauge dependence. The colour group Casimirs 
are defined by 



rpdrpa 



Cf 



jacd jbcd 



CaS' 



ab 



Tr r"r^ 



Tpd' 



ab 



(2.9) 
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where are the colour group generators with structure constants f"'^'^ and 1 < a < Na with A^^ 
the dimension of the adjoint representation. The coupling constant appearing in the covariant 
derivative is g but appears in (12. 8p in the combination a = j (IGvr^). The absence of a correction 

io'i^'?' [a] to this order is a reflection of the all orders result where there is no renormalization 

III \ ) 

of dWi. This is because it corresponds to the vector current which is a physical operator and is 
not renormalized reflecting the Slavnov- Taylor identity of the underlying gauge theory. See, for 
example, [35] for background. Moreover, this property is retained in all schemes and we must 
ensure that the RI'/SMOM scheme renormalization respects this general feature. 

In order to define the RI'/SMOM renormalization scheme for level Wi we compute the 
Green's function illustrated in Figure 1 where the external momenta p and q are independent. 
However, we will focus totally on its structure at the symmetric subtraction point, [27], which 
is defined as 

/ = = + = - (2.10) 
where \x is the common mass scale. These relations imply 

Vq = \lJi^ . (2.11) 

In order to determine the renormalization constant in the MS or RI' schemes one can proceed 
with only one external momentum. For instance, the choice p = — q corresponds to a zero 
momentum insertion before setting p'^ to be a specific scale involving /i^. Though in fact to de- 
termine the full mixing matrix in the MS scheme requires accessing the total derivative operators 
which is achieved by nullifiying the momentum of one of the external quark legs, [32] . As the 
operators we consider have free Lorentz indices then the Green's function has to be decomposed 
into a set of Lorentz invariant amplitudes. The basis for W2 has already been constructed in 
[31] and we briefly recall the construction as well as the notation. First, we write the Green's 
function as 



n,: 



(^(p)Oi......„X-p-9)^(«))| . , . = T.'Pw,....,.SP^^)^(kM'l)- (2-12) 

^ ^ ^ k=l 

Here S^'^ {p, q) are the scalar amplitudes for each of the two operators O* whereas P^*^^ {p, q) 

are the Lorentz tensors of the basis whose Lorentz indices reflect the symmetries of the operator 
which has been inserted. For level W2 there are ten independent Lorentz tensors and we will 
detail these later. To evaluate the scalar amplitudes we use the method of projection where we 
apply different linear combinations of the basis tensors to the Green's function in such a way 
that each amplitude is isolated in turn. More speciflcally we have, [32], 

^fk)iP,q) = MhVl,{^-'^-{p,q) ((v(p)0;.,...,„,(-p-g)V5(5)))|^,^^,_^, (2.13) 

where A^^; is a matrix of rational polynomials in the spacetime dimension d. It is constructed 
by first determining the matrix Mj^i which is defined by the contraction of the tensors of the 
basis with each of the other tensors via 

j^ii = n^)......„,(^'>'?)^ar'"'"(P'«)L.=,.=_,. • (2.14) 

Then is the inverse of AA^^. 

More specifically, the basis of Lorentz tensors into which we decompose the Green's functions 
with the operator insertions is. 



'P^f^.uiP, q) = l^iPu + luP^, - ^hi^u , V^)^^{p, q) = -ff,qu + Juqt, - ^i^t^v 
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V. 



V. 



W2 

W2 



V, 



W2 



1 1 

■j^Pi^Qiy + -^Qi^Pf 



V. 



V. 



W2 
W2 



-J'^^lu 



{p,q) = i) 

(.p,q) = i 



■ 1 1 1 
1 1 

-J^pf^p'^ + ^^^^ 



vZl..,{p.q) 



T)W2 

' (10)/ii/ 



(p, q) 



1 
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^{3)fj.pqPi' + ^i3)upqP^l 
r(3) fipqqu + r(3) ,ypqq^ 



(2.15) 



As there are ten elements in the basis then in order to save space for presenting the projection 
matrix, we have partitioned the 10 x 10 matrix into four sub-matrices. Defining 



M 



W2 



1 



W8{d - 2) 



(2.16) 



then we have 



/ 18 



MY,' 



48 
24 
12 



9 

18 
24 
30 
24 



24 
12 

32(d + 1; 
16(d + 1; 
8{d + 4) 

/ 24 12 



48 
24 

64(d + 1) 
32(d + 1) 
16(d + 4) 

30 
24 

) 16(d + 4) 
) 20(d + l) 
16(d + 1) 



24 

30 
32(d + 1 
8{5d - 1 
8{4d + 1 

24 
48 

8(d + 10) 
16(d + 4) 
52{d + l) 



mY^ 



30 
24 





24 
48 





mZ' 



( 32(2(i 
8(4d + l) 
16(d + 4) 


V 



32(d + 1) 
16(d + 4) 
8((i + 10) 





8(4d + 1 
8(5d - 1 
32(d + 1 





1) 



16(d + 1) 
20((i + 1) 
16((i + 4) 



16(d + 4) 
32(d + 1) 
64(d + 1) 





12 \ 

24 

16(d + 4) 
8(4d + 1) 
32(2d-l) j 

\ 




/ 

8(d + 4) 
16((i + 1; 
32((i + 1; 









-24 
-12 







-12 
-24 / 



(2.17) 



We note that here the label W2 refers to the level and the same projection is used for dW2- In 
defining this basis we have used the set of generalized 7- matrices denoted by F^^j"'^" which are 
totally antisymmetric in the Lorentz indices and defined by, [361 EH EH] j 
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•7' 



(2.18) 



where the notation includes the overall factor of 1/n!. For simplicity we will retain the more 
natural notation for a single 7-matrix rather than use the clumsy Fj^-^^ . General properties have 
already been discussed in [391 HO] but we note that one of particular use is 



tr 



m) (n) 



OC 



5n 



..fl,-al/l...l/n 



(2.19) 
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where I^^ 



../lmI^l..Mr 



is the unit matrix in the space of these generahzed 7-matrices. This leads to 



the natural partition of Ai^'^ into two submatrices. The use of F^^j"^" is especially appropriate 
as we will be using dimensional regularization and these objects form the complete basis for 
spinor space in d-dimensional spacetimes. Moreover, they allow us to see that the set (|2.15p is 
complete. If there were more than two independent momenta in this problem then one would 
have to include matrices with n > 4. One final comment on the basis choice and that is that the 
basis is not unique. One could choose different linear combinations of the tensors we use here 
but the overall structure of the Green's function would be unaltered. Though we have chosen 
to retain a degree of symmetry if one interchanges p and q. For the total derivative operator, 
dW2, which is symmetric under this then the amplitudes will also respect this property as will 
be evident in the explicit expressions. 

Having dealt with the general structure of {ip{p)0\ „ {—p — Q)i^{Q)) 00 ^ now 

need to outline the procedure to renormalize our operators in the RI'/SMOM scheme. There 
are a variety of ways of doing this. If for the moment we denote by a particular amplitude or 
combination of amplitudes which contains all the divergent parts of the Green's function then 
we will use the renormalization condition 



lim 



1 



(2.20) 



where d = 4 — 2e in dimensional regularization. This is in keeping with the overall regularization 
invariant scheme ethos that there should be no 0{a) finite parts for this particular amplitude 

combination 0, [271 [28l [29] . The renormalizat ion constant is the one which removes 

the divergences and leads to the anomalous dimension of the operator. However, as the operator 
has been inserted into a quark 2-point function then the wave function renormalization constant 
of these external quarks has to be considered which is reflected in . The annotation here is 
RI' as the 2-point functions of the theory are renormalized according to the original procedure 
of [211 [22| [23] . In the RI' scheme the 3-point and higher functions are renormalized according to 
the MS scheme so that their finite parts after renormalization are not unity. In [23j the Landau 
gauge three and four loop RI and RI' scheme quark field and mass anomalous dimensions were 
recorded. The anomalous dimensions of all the fields in an arbitrary linear covariant gauge were 
given at three loops in [24]. In the latter article the renormalization of the gauge parameter 
was also determined and the relation between the parameter definition in different schemes was 
established. For instance, if we denote the scheme in which the variable is defined in by a 
subscript then 

«RI' = "MS + O (a|js) (2.21) 

and, 



"RI' 



+ ((l8a|^ - 18a|jg + 190a|^ - 576C(3)aj^ + 463aj^ + 864C(3) - 7143) 
+ (- 320a|jg - 320aj^ + 2304C(3) + 4248) CaTfNj 

O -1 



+ 



(- 4608C(3) + 5280) CpTrNf) ^ + O (^a^^ 



(2.22) 



These relations are important as we will be providing results in arbitrary linear covariant gauge 
in each scheme separately. Therefore the conversion of parameters between schemes has to 
be taken into account when comparing results and, moreover, when constructing conversion 
functions for the operators. 
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One reason we provide the amplitudes in both the MS and RI'/SMOM schemes is that it 
turns out that there is no definitive scheme definition for the latter. This is because in discussing 
the general scheme definition we alluded to channel which represents some combination of the 
scalar amplitudes. This choice is not unique as indicated in [31]. For instance, one could 
extract the combination given by multiplying the Green's function by the Born term which was 
the approach of [271 ES] for the tensor current. However, one could equally ensure that that 
amplitude with the poles in e has no 0(a) terms after renormalization. Indeed this statement 
does depend on how one defined the basis tensors and therefore one has a range of choices as to 
how to define RI'/SMOM. Though as an aside the renormalization of dW2 has to be carried out 
in such a way that it is consistent with the Slavnov- Taylor identity, |31j . since this operator is 
related to the divergence of the vector current which is conserved and a physical operator with 
zero anomalous dimension to all orders in perturbation theory in all renormalization schemes. 
The version of the RI'/SMOM scheme which we will use here is a direct extension of the one 
loop version given in [31]. There the coefficients of the channels 1 and 2 were used to define 
the W2 renormalization constant. Two channels were used due to the fact that we have two 
renormalization constants to fix for the first row of the mixing matrix of ()2.4I) . As with the 
three loop MS renormalization of [32] the counterterms are entwined with each other and so 
one has to solve a set of linear equations to fix the explicit values for each renormalization 
constant. Therefore, all the amplitudes which we record for RI'/SMOM are determined via 
that condition. Though we emphasise that this choice is not unique in order to define the 
renormalization constants. One could have, for instance, multiplied the Green's function by 
the Born term. Alternatively one could use the same combination of amplitudes as was used 
for dW2, [31j, for consistency with the Slavnov- Taylor identity of that operator. Though that 
would not be sufficient on its own as one would require a second independent projection to 
solve for both elements of (12. 4p . However, in order to assist with lattice computations in the 
situation where a different amplitude combination might be made to renormalize the operator, 
we provide the MS amplitudes as it is the canonical reference scheme as well as to facilitate 
making different variations on the scheme definition. This is because one can then readily 
convert from any scheme to MS for comparison. Indeed some lattice groups carry out their 
renormalization according to their own prescription before converting their results to MS prior 
to performing the actual matching to the continuum results. 

Finally, having described the theoretical background to the problem we note how it was 
implemented in practical terms. The main tool used to handle the algebra of Lorentz indices, 
projection matrices and evaluation of the underlying Feynman graphs was the symbolic manip- 
ulation language Form, [41j . The actual Feynman diagrams themselves were constructed in 
electronic format by using the Qgraf package, [32]. These were then converted into the nota- 
tion used for the Form computations where the colour and Lorentz indices were added. At one 
loop there were 3 graphs and at two loops there were 37 diagrams. The algorithm to evaluate 
the integrals comprising each Feynman graph was to first write them as scalar integrals. By 
this we mean in a form which is the starting point of the Laporta algorithm. Each integral is 
rewritten in terms of the denominator propagators as far as possible. For those cases where this 
is not possible the numerator tensor structure is written in terms of irreducible propagators. In 
this form one can apply the Laporta method, j43j , where a redundant set of linear equations is 
established between all the integrals by using integration by parts and Lorentz identities. These 
relations are then solved to express each required integral in terms of a small set of scalar master 
integrals. These have been evaluated directly and given in a set of articles, [441 [55l H6l I47j . A 
summary has been provided in ^29j. Central to the construction of the Laporta algorithm used 
for the current problem was the use of the Reduze package, [35], which is based on the GiNAC 
system, |49) . which is written in C-|— 1-. For the two loop Feynman integrals we needed to eval- 
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uate, it turns out that there are only two basic momentum topologies for which one needs to 
construct the reduction of integrals. These are the ladder and its associated non-planar ladder. 
The momentum routing of all the Feynman graphs could be mapped to these two topologies 
and therefore, aside from the elementary one loop case, Reduze was only required to build two 
sets of integration by parts and Lorentz identity relations. The machinery derived for this set 
of operators was also used in constructing the two loop amplitudes for the quark currents, [30] , 
There the anomalous dimensions for the scalar and tensor currents reproduced the known results 
of \27\ [28l [29] and therefore this provides us with a strong check on the routines which have 
been built and used for level W2 here. Finally, in order to carry out the renormalization prior to 
deducing the values of the amplitudes, we follow the method of [50] for automatic calculations. 
We perform the computation for bare parameters and operators and then introduce the coun- 
terterms by rescaling all bare quantities with their renormalized equivalences. For the operators 
themselves this ethos also applies in that the bare operators are replaced by their renormalized 
counterparts defined by the structure of the mixing matrix, ()2.4p . Then the counterterms for 
the operators are fixed with respect to whichever is the scheme prescription after the known 
counterterms for the coupling constant and gauge parameter have been included from previous 
calculations. 



3 Amplitudes. 

This section is devoted to recording the amplitudes for the two second moment flavour non- 
singlet Wilson operators, W2 and dW2, inserted in the quark 2-point function. As there are 
quite a large number of amplitudes and two operators to provide the full set of expressions for 
in each scheme, we present the results of the finite parts in a set of Tables. We follow the same 
notation as [30j in recording the sets of numbers which appear as coefficients of both the group 
Casimirs as well as a set of basis numbers for each amplitude. The latter derive from the explicit 
values of the basic two loop scalar master Feynman diagrams which were summarized in |29) . 
More concretely our amplitudes are written in general as 

+ [T.-iT'-'n'^)cFCAa^ + (E^S'r«i^^^)c^V + O(a^). (3.1) 

Here the On^ are the basis of numbers which for simplicity also includes the gauge parameter. 

The label / here indicates both the loop order, as the first number in the two loop case, whilst 

(I) 

the second number at two loops relates to a specific colour group Casimir. The coefficients c^^-^'^ 
are the actual rational numbers which appear in the appropriate piece of the amplitude as is 
evident from the Tables themselves. Clearly these reduce the space needed to display the explicit 
results and allows for easier comparison of the structure of the amplitudes between schemes. For 
the numbers in the basis we use the notation of [29]. For instance, il^{z) is the derivative of the 
logarithm of the Euler Gamma function and 

(3.2) 

where Lin{z) is the polylogarithm function. The quantity S is a combination of various harmonic 
poly logarithms, [SSI HZ]; 

s = n'i^ + n^S (3.3) 
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and this combination always appears. 

Whilst the Table£] reflect the size of the computation undertaken for more practical purposes 
it is appropriate to provide the explicit numerical evaluation of the results. In order to achieve 
this we note that the explicit numerical values of the various numbers in the basis are taken to 
be 



C(3) = 1.20205690 , E = 6.34517334 , V' (\ 



S3 



,1838167 , S2 



0.32948320 , S3 



vr 



vr 



= 0.32225882 , S2 
0.19259341 . 



10.09559713 

= 0.22459602 



(3.4) 



Therefore, we can record the numerical values in both schemes. We have chosen to do this for 
the case of the SU{3) colour group. For the MS scheme we have 



MS 



MS 



+ 



[1.7499534 + 0.4444444a] a 

37.3849283 + 2.4296818a + 1.6608855a2 - 5.0243682iV/] + O(a^) 
1 + [3.3748835 - 0.1387491a] a 
43.5097605 - 0.7932193a - 0.1511137a2 - 4.7881095iV/l + O(a^) 



_ = [1.2531175 + 0.2037969a] a 



+ 



15.7237696 + 3.9484943a + 0.7217806a^ - 0.9599542A(^j + 0{a' 



_ = [1.6541967 + 0.4075938a] a 



MS 



+ 



16.8558938 + 4.8339947a + 1.4319988a2 - 1.7451297A(/ 



[2.6233015 + 1.7040764a] a 



+ 



60.6554860 + 15.1801309a + 6.1487315a2 - 4.8269741iVj 



_ = [0.9322541 + 0.6850920a] a 



+ 



28.1056849 + 9.0379699a + 2.4124553a2 - 0.90U209Nf 



c? + 0{c?) 



_ = [1.5956635 + 1.0926858a] a 



MS 



MS 



MS 



+ 



+ 



46.3221825 + 11.8708374a + 3.9440029a2 - 2.4490197iV/j 0? + 0{o? 
[1.6910474 + 0.4075938a] a 

21.8712121 + 5.2430694a + 1.4690359a2 - 1.6999495A^/1 + Oic?) 
0.4444444a 

8.8385249 + 1.3267080a + 0.0370370a2 - 0.4428779iV/ 

- 1.6390287a 

- 30.9488446 - 0.9782418a + 0.1365857a^ - 2.5665832iV/ 
+ 0(a3) . 



c? + 0(a= 



(3.5) 



* Attached to this article is an electronic file where all the expressions presented in the Tables are available 
a useable format. 
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Those for dW2 are 



MS (2) ^^'^^ 



MS 



- 1 + [1.6249301 - 0.5831936a] a 

+ [6.1248321 - 3.2229010a - 1.81199920^ + 0.2362586iV/ 



MS 



^dW2 
^(8) 



{p,q) 



MS 



[2.9441649 + 0.6113907a] a 



+ 



37.5949817 + 9.1915636a + 2.1908165a2 - 2.6599037iV/ + O(a^) 



= s 

MS 



dW2 
(7) 



{p,q) 



MS 



[3.2498602 + 1.5002795a] a 



+ 



63.1780763 + 16.7048322a + 5.37600170^ - 4.1941494iV/ + 0(a 



^(5) 



MS 



^dW2 



iP,Q) 



MS 



[3.5555556 + 2.3891684a] a 



+ 



^.7611709 + 24.2181008a + 8.56118690^ - 5.7283951iV/ + 0(a 



MS 



^(10) 



MS 



- 2.0834731a 

- 39.7873696 + 0.3484662a + 0.17362280^ - 3.0094611iV/ 



(3.6) 



The exphcit results from which these numerical values were derived are recorded in Tables 1 
to 12. However, the equivalence between amplitudes indicated above for dW2 are exact at 
two loops which is why we have not included parallel columns with the same values. Indeed 
these equivalences are a minor check on our computation as the tensor basis is clearly reflection 
symmetric when the operator insertion has the same property as is the case for dW2- Another 
check resides in the fact that some of the amplitudes for dW2 are proportional to the two loop 
amplitudes for the vector current of [3^. For instance, 



^(1)^ 



MS 



MS 



MS 



^(2)^ 



_ = is 

MS 



dW2 
(3) 



(p, q) 



MS 



MS 



MS 



MS 



(3.7) 



where V indicates the vector current in the notation of [3T1 132] . Though it should be stressed 
that the associated amplitude basis tensors of each of the operators are not the same. This is 
trivial to see because of the mismatch in the number of Lorentz indices on each of the operators. 
Indeed this is why channel 4 of dW2 does not feature in p.7p as this index imbalance means 
that there are not the same number of amplitudes for each operator. 

To illustrate that the Slavnov- Taylor identity is satisfied by the operator dW2 we have to 
project out that part of the Green's function which contains the divergence of the vector operator 
and was discussed in [30\ [31] . This is achieved by simply contracting the two free indices of 
the dW2 operator. Prom the explicit forms of the amplitudes in the Tables the combination 
proportional to ^ is 



+ 



[d-2] 
d 

[d + 2] 
2d 



^dW2 



(p, q) 



MS 



MS 



+ 



^(2) 

[d-4] 
Ad 



MS 

^dW2 
^(5) 



+ 



[d-4] 
Ad 



^t^'{P,q) 



MS 



(p, q) 



MS 



12 



- + -Cpaa 




4 lb 



(3.8) 



Clearly the right hand side is proportional to the finite part of the quark 2-point functions after 
renormalization in the MS scheme. A similar combination gives the part proportional to ^ with 
the same expression to two loops. We emphasise that the MS renormalization constant used to 
renormalize dW2 is the same as was used in [32j to construct the operator correlation functions. 

For the RI'/SMOM scheme we have to be careful to define the renormalization of dW2 so 
that the Slavnov- Taylor identity is satisfied. In other words the renormalization constant for 
dW2 is already determined since the vector current is a physical operator and therefore its 
renormalization constant is unity in all schemes. However, for W2 itself one has a large degree 
of freedom to define the RI'/SMOM scheme renormalization constant. This was discussed in 
|30l [3T| and we follow the prescription used there. This was to fix the finite parts of the other 
renormalization constants of the W2 level mixing matrix so that for those channels containing 
the poles in e there was no 0{a) corrections. Of course this is not the unique way of defining 
these renormalization constants. One could, for example, take some sort of projection of the 
Green's function and require that there is no 0{a) correction for that particular combination. 
The fact that there is a degree of ambiguity is one of the reasons why we have provided the 
MS results. Given that we are extending [31] to two loops we follow that prescription here and 
record those results. The complete results are in Tables 13 to 23 and the SU{3) numerical results 
for the operator W2 itself are 



[1 



,2531175 + 0.2037969a] a 
16.9936409 + 4.6449275a + 0.81800470^ + 0.15284770^ 

- (0.9599542 + 0.2264410a)iV/] + O(a^) 
,6541967 + 0.4075938a] a 

19.6462784 + 6.6534412a + 1.50119810^ + 0.3056953o^ 

- (1.7451297 + 0.4528820a)iV/] + O(a^) 
6233015 + 1.7040764o] a 

67.8774218 + 32.7286131o + 8.16392170^ + 1.2780573a^ 

- (4.8269741 + 1.8934182a)iV/] + O(a^) 
,9322541 + 0.6850920o] a 

26.6612784 + 12.6105871o + 2.58767010^ + 0.51381900^ 

- (0.9014209 + 0.7612133a)iV/] + O(a^) 
,5956635 + 1.0926858o] a 

48.8125936 + 22.7212025o + 5.2502695o^ + 0.8195143o^ 

- (2.4490197 + 1.2140953o)iV/]a2 + O(a^) 
,6910474 + 0.4075938o] a 

25.3854030 + 8.7651597o + 1.93329690^ + 0.3056953o^ 



+ 



[1 



+ 



+ 



+ 



[1 



+ 



[1 



+ 
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- (1.6999495 + 0.4528820a)iV/] + 0{a^ 
O.UUUia 

7.4702494 + 0.5365824a + 0.0370370a2 

1.6390287a 



0.4428779iV/ + O(a^) 



- 35.7026184 - 1.4031864a + 0.1365857a2 - 2.5665832iV/ 
+ Oia^) . 

Those for the total derivative operator are 



\2) 
1 + 



[1.6249301 + 0.7501398a] a 



+ 



■^8W2 
^(3) 



31.5890382 + 12.2494724a + 2.8130241a^ + 0.5626048a^ 
- (2.0970747 + 0.8334886a)Af/] + 0{a^) 

{p,q) 



^(8) 

[2.9441649 + 0.6113907a] a 



+ 



^(4) 



(p, q) 



37.5949817 + 10.2080849a + 2.2927149a2 + 0.4585430a^ 
- (2.6599037 + 0.6793229a)iV/] + O(a^) 

{p,q) 



^(7) 

[3.2498602 + 1.5002795a] a 



+ 



^dW2 
^(5) 



(p, q) 



63.1780763 + 24.4989449a + 5.6260483a2 + 1.1252097a3 
- (4.1941494 + 1.6669773a)iV/]a2 + O(a^) 

{p,q) 



^(6) 

[3.5555556 + 2.3891684a] a 



+ 



88.7611709 + 38.7898049a + 8.9593816a2 + 1.7918763a3 
- (5.7283951 + 2.6546316a)iV/] + O(a^) 



^dW2 
^(9) 



(p, q) 



- 2.0834731a 

- 39.7873696 - 2.4294979a + 0.1736228a^ - 3.0094611iV/ 
+ 0(a3) . 



(3.9) 



(3.10) 



The same remarks made for the MS case concerning the equivalences indicated above between 
amplitudes for dW2 and their proportionality with those of the vector operator of [30] also apply 
to the RI'/SMOM scheme results. So the relations of (j3.7p are true when the MS indication is 
dropped. Equally we note that the Slavnov- Taylor identity is also satisfied in the RI'/SMOM 
scheme. The parallel computation to (|3.8p is 



+ 



[d-2], 

d ■ 
[d + 2] 



{p, q) 



^dW2 
^(2) 



{p,q) + 



[d - 4] ^dW2 



(p, q) 



^'^HP,q) + ^^^'^Hp,q) 



I + 0{a') 



(3.11) 



2d W ^^'^^ ' Ad (5) 
to the order we have computed to. As the quark wave function renormalization constant is 
in the RI' scheme there are no 0(a) corrections which indicates consistency with the Slavnov- 
Taylor identity. As a final observation on the two loop values of the amplitudes we note that 
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the relation noted in |3lj at one loop 

is valid to two loops for all a and for both the MS and RI'/SMOM schemes. 



(3.12) 



4 Anomalous dimensions. 



In order to record these amplitudes we have determined the renormalization constants which 
were fixed by the RI'/SMOM scheme choice defined in [31]. These are encoded in the mixing 
matrix of anomalous dimensions which extends (12.81) and is 



RI'/SMOM 



-Cpa 



+ 



^lOSa^ + 324a - 924)^^ 



[12o? + 216q; - 616)7r^ 



81q^ - 243a + 16866 



+ 



Ca - 2016C^ 
Cfc? 



336^'(i) - 2247r^ - 5976 TpNf 



486 



+ 0{c? 



7i^'(a,a) 



RI'/SMOM 



-Cpa 
3 



+ 



264V'' 



3> 



Ca 



7^' (a, a) 



RI'/SMOM 



0(a' 



1767r^ - 81a^ - 243a - 6651 

'^i^a - 288 + 648a - 967r^ + 1927r^a 



144V''(i) - 288V'\3 



2„2 



647r^ + 2340 - 96^'(i) TpNf 



486 



+ O(a^) 



Cf 



(4.1) 



As a check on the derivation of these we ensured that the two loop MS matrix of anomalous 
dimensions emerged from the same FORM programme. They agreed with the original computa- 
tion of ^32j which provided the matrix at three loops. However, that computation was performed 
with the use of the Form version of the Mincer algorithm, [5HI52|. There one could deduce the 
off-diagonal matrix element by choosing to route the single external momentum out through the 
operator insertion itself in contrast to the current computation where there are two independent 
momenta. As dW2 is in essence the vector operator then the 22 element is equivalent to the 
vector operator anomalous dimension. 

The RI'/SMOM scheme anomalous dimensions can be computed by a second method in- 
volving conversion functions. For an introduction to these see [35] • However, for level W2 this is 
not as straightforward since one is not dealing with a multiplicatively renormalizable operator. 
Instead there is a matrix of renormalization constants and therefore the concept of a conversion 
function translates into a matrix of conversion functions. Formally, for W2 this is 



a) 



7W2 
ifc.RI'/SMOM 



7W2 

fcj.MS 



(4.2) 



The explicit forms of the elements of this matrix are straightforward to deduce and are 

r^W2 

„W2( ^ ^ll,RI'/SMOM W2( ^ 



71^2 
22,RI'/SMOM 



y^W2 

ll,MS 



2;W2 

22,MS 



15 



^Wif \ 12,Rl7SMOM 11,RI7SM0M i2,MS on 

Ci2'(«>«) = -^w^ yW2yW,_ (4-3) 

22,MS 11,MS 22,MS 

where (a, a) = from the upper triangular nature of the underlying renormalization ma- 
trices. Indeed this means that the diagonal conversion functions are what one might naively ex- 
pect. From this one can derive the relation between the anomalous dimensions in both schemes. 
Specifically we have 
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nil'/SMOM («RI'' «RI') - ^u'ms ("ms) " 1^ («Ms) '^n' («MS' "Ms) 



«MS^^^ («MS' «Ms) 1" ("MS' "Ms) (^.4) 



W2 { \ 

^12,RI7SM0M i"RI''"RI'i 



and 



7. 



W2 

22,Rl7SMOM 



iORI''«RI' 



7. 



W2 

12 



( \ ^W2 ( A 

MS V^MsJ '-^11 v"ms'"msJ 

^ («Ms) a|=^it (- 



MS'^MSJ 



•^MS'^MS 



("■MS'^Ms) f)r, '^12^ ( 

VI/2 Z' \ ,^14^2 ( A 

'^ll,MS V'^MSJ '-^IS 1,«MS'"msJ 

^22^8 ("ms) '^12' ("mS' "ms) 
+ (a^S' "Ms) ^ ("Ms) "^n' ("MS' "Ms) 

+Clt( 



+ 



"mS'^MsJ "MS^a 



da 



'MS 



ms'"ms; 



("ms'"ms) 
(^ 



'ms'"ms; 



(4.5) 



W2 

'22 



HS ("Ms) - P («Ms) ^-2^' («MS ' " 

Ms^o^^ ( 



\J_ 

"MS'"MSJ Q(y. 



— lnC2^^(aj^,a^) (4.6) 
MS 



where to avoid confusion we have labelled the scheme the variables are in explicitly. We have 
used the designation RI' for the variables on the left side of the equations as we use the definitions 
of [23] which were derived from renormalizing the 2-point functions of all the fields. With these 
definitions we have computed the conversion function matrix explicitly to two loops and fincO 



1 + 



(36a - 42)^'(i) + (28 - 24a)7r2 - 27a + 459 



+ 



(5184a^ - 12096a - 4608)(^'(i))2 + (16128a - 6912a^ + 6144)^'(i)7r' 



+ (328536a - 13608a^ - 613656)V''(5) - (5022 + 1944a)V'"'(;^ ' 
+ (1259712a - 5248800)s2(f ) + (10497600 - 2519424a)s2(|) 
+ (8748000 - 2099520a)s3(f ) + (1679616a - 6998400)s3(f ) 



^These together with the anomalous dimensions are also included in the attached electronic file. 
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+ (23040^ - 192a + 11344)7r^ + (9072a^ - 219024a + 409104)7r^ 
+ 7290a^ - 90882a + 107568 + (34992a - 40824)S 

ln2(3)7r 



+ (104976a + 559872)C(3) + (8748a - 36450)- 



+ (437400 - 104976a) i^^^^ + (39150 - 9396a) ^ 

V 3 v 3 



V3 

Cf 



+ [5832(V''(|))^ - 7776V''(^)7r2 + (11664a2 - 39366a + 99468)V''(^) 

+ (5103 - 972a)V'"'(|) + (2519424 - 629856a)s2(f ) 

+ (1259712a - 5038848)s2(f ) + (1049760a - 4199040)s3(|) 

+ (3359232 - 839808a)s3(5) + (2592a - 11016)7r^ 

- (7776a2 - 26244a + 66312)7r2 - 8748a2 - 34992a + 1759644 

+ (17496a - 34992)S + (26244a - 734832)C(3) 

+ (17496 - 4374a)i^^^^^ + (52488a - 209952)-''^^^'' 



+ (4698a - 18792) — 
v3 



^/3 V2, 
Ca 



+ 



46656^^' (I) - 3110471^ - 642168] TpN^ + O{o?) (4.7) 



C}^2(a, a) = 24V''(3) - 167r^ - 54a - 297 



; Cpa 
162" 



+ 



13824a + {\)f - (19440a^ + 18432a + A&M%)ii'{\)T{^ 

- (342144a + 6912a2 - 824904)V^'(i) + 18792V'"'(|) 

+ (10917504 - 1679616a)s2(f ) + (3359232 - 21835008a)s2 (f ) 
+ (2799360 - 18195840a)s3(f ) + (14556672 - 2239488a)s3(f ) 
+ (6144a - 34496)7r^ + (12960a2 + 228096a - 549936)7r2 
+ 37908a^ + 76788a - 188892 + 46656S 

- (139968a + 1609632)C(3) + (75816 - 11664a) ^ 

v3 

+ (139968a - 909792)^^ + (12528a -81432)^ Cp 

v3 v 3 



+ 



31104V''(|)7r2 - 23328(V/(^))^ + (107892a - 5832a2 - 122472)V''(^ 
+ (972a - 9720)V''"(i) + (944784a - 4513968)s2(|) 
+ (9027936 - 1889568a)s2(f ) + (7523280 - 1574640a)s3(f ) 
+ (1259712a - 6018624)s3(f ) + (15552 - 2592a)7r'^ 
+ (3888a^ - 71928a + 81648)7r2 - 21870a^ 
+ 34992S + (1277208 - 52488a)C(3) 
ln2(3)7r 



113724a - 2488482 



+ (6561a -31347) - 



V3 



+ (376164 - 78732a) 



ln(3)7r 



+ (33669 - 7047a) 



TT 

71 



Ca 
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+ 

and 



Cwa? 

293767r^ -44064V'' (i) +946080 TpNf f ^ + O(a^) (4i 

104976 



Cf2'{a,a) = 1 + 0{a^) . (4.9) 
For SU{2>) the numerical values are 

CYi{a,a) = 1 + [1.6390287a + 5.1248369]a 

+ [6.51081510^ + 23.1781730a + 132.6228486 - 12.1458110A('/-]a2 + 0{a^) 
C7j^2(a,a) = - [0.4444444a + 1.7499534]a 

- [2.1301455a^ + 7.2772405a + 49.1967063 - 5.0243682Af/]a^ + O(a^) 
Cf^^a^a) = 1 + O(a^) . (4.10) 

As the operator dW2 and the vector current are in effect equivalent then both their anomalous 
dimensions and conversion functions are the same. For the latter case this is the 22 element of 
both matrices. 

With these values we can now make a comparison with similar functions in the RI' scheme 
for the level W2- However, in order to do this we need to define the conversion functions 
from the renormalization constants. This is not as straightforward for RI' since the momentum 
configuration defining the scheme does not access the off diagonal elements of the full operator 
mixing matrix. Given that we have chosen to define our operator basis in such a way that this 
matrix is triangular it is possible to define and compare the diagonal elements of the conversion 
function matrix. As noted for RI' only the diagonal elements of the renormalization constant 
matrix exist and so we define the two conversion functions of interest as 



= , C'2t(a,a) = -1^. (4.11) 

ll,MS 22,MS 



The explicit forms were given in I25j and in the Landau gauge we have the SU{2>) numerical 
values 

CYi{a,Q) = 1 + 4.5925926a + [119.8268158 - 10.9794239iV/]a^ + O(a^) 
C2^^(a,0) = 1 + 0{a^) . (4.12) 

Clearly the numerical values of the one and two loop corrections for the RI' scheme are each 
about 10% smaller than those for the RI'/SMOM scheme for Nf = 2 and 3. However, it is 
not entirely clear whether this is a proper comparison in the sense that one is not necessarily 
comparing with respect to the same tensor basis. Equally the mixing matrices are not truly the 
same as there can be no off-diagonal element for RI'. Indeed it is probable that the convergence 
of the RI'/SMOM scheme result could be improved by a different choice of basis tensors or 
exploit the freedom in how one defines the RI'/SMOM scheme for this operator. That aside one 
does at least have the explicit forms of the amplitudes at the symmetric subtraction point at 
two loops in order to assist with lattice matching. 

One benefit of the explicit forms of the conversion functions is that we can deduce the three 
loop RI'/SMOM scheme anomalous dimensions in the Landau gauge using the three loop MS 
expressions of ()2.8p . We have 



7iT^(a,0) 



RI'/SMOM 3 
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+ 



+ 



SOStt^ 



462V''(i) + 84331 Ca 



+ 



IOOSCf 
Cpa^ 



168V;'(i) - 1127r2 - 2988] TpNf] + 0{a') 

64152(V''(|))^ - 85536^/;'(i)7r2 + 862812^'(i) 

+ 56133V'"'(|) + 27713664s2(f ) - 55427328s2(f ] 

- 46189440s3(|) + 36951552s3(f ) - 1211767r^ - 

- 384912S - 7243344C(3) + 25273296 + 192456 ^ y 

V3 



5752087r2 

2/ 



- 2309472 



+ 



+ 



+ 



+ 



lii(3)7r 



- 206712^ 

V3 



119328V''(5)7r2 - 89496(^'(|))^ - 5901984V''(5) 

- 55242V''"(i) - 57736800s2(|) + 115473600s2(f ) 
+ 96228000s3(f ) - 76982400s3(f ) + 1075367r^ 

449064S + 3639168C(3) - 4833270 

.3 



+ 39346567r^ 

,2 



- 400950^^^ + 4811400i^ + 430650^ 
v3 v3 v3 



CaCf 



31104(^'(i))2 



23328V'' (i)vr2 + 251424V''(|) 



20412V' 



/" ( 1 



10077696s2(f ) + 20155392s2(f ) 
+ 16796160.S3(|) - 13436928s3(f ) + 440647r'' 
- 1676167r2 + 139968E + 1259712C(3) - 16603056 



69984^^^5}^ + 839808^ + 75168 " 



V3 



V3 



V3 



CATpNf 



32544(V''(i))^ - 43392V'' (i)7r2 + 2227824Vy(i) 

+ 20088V'"'(i) + 20995200s2(f ) - 41990400s2(f ; 

- 34992000s3(|) + 27993600^3 (f) - 391047r^ 

- 14852167r^ + 163296S - 559872C(3) - 742608 



+ 145800 



ln2(3)7r 



1749600 



lii(3)7r 



V3 ^ 
1244167r2 - 186624V'' (I) + 2423520 



- 156600^ 
V3 



,2 .j2 



+ [1679616C(3) - 90720] Cp 



Cpa^ 
39366 



CFTpNf 



(4.13) 



Upa 

3 



+ 



+ 



+ 



264V''(i) - 1767r2 - 665ll Ca + [l44V''(i) - 967r2 



288 



Cf 



647r^ - 96V' (I) + 2340 TpNf 



486 



342144V'' (i)7r^ - 256608(V''(^))^ - 1082808V''(^) 

- 106920Vy"(i) - 49653648s2(f ) + 99307296s2(f ) 

+ 82756080s3(f ) - 66204864^3 (|) + 1710727r'^ + 7218727r2 
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and 



+ 384912S + 12369672C(3) - 37423134 - 344817 



l n^(3)7r 



+ 4137804i^^^^5^ + 370359^ 



Cl 



+ 



+ 



[477504(^'(i))2 - 636672^' (i)7r^ + 7978176V''(5) 

+ 204768V'"'(|) + 117993024s2(f ) - 235986048s2(f ) 

- 196655040s3(f ) + 157324032s3(f ) - 3338247r^ 

- 53187847r^ + 653184S - 11897280C(3) + 367416 

+ 819396^^^^ - 9832752^^ - 880092^1 CaCf 
V3 V3 VsJ 

[93312(V''(i))^ - 124416V'' (Dtt^ - 150336V''(|) 

+ 38880V''"(5) + 18055872s2(f ) - 36111744s2(|) 

- 30093120s3(f ) + 24074496s3(|) - 622087r'' 

+ 1002247r^ - 139968S - 1749600C(3) + 24312312 



+ 125388 



ln2(3)7r 



+ 



x/3 

208896V'' (i)vr2 - 
- 75168V'"'(i) 



1504656 



ln(3)7r 



134676-= 

\/3 



CATpNf 



156672(V''(i))^ - 3297024V''(i) 

- 43670016s2(f ) + 87340032s2(f ) 
+ 72783360s3(f) - 58226688.S3 (f ) + 1308167r'' 
+ 21980167r2 - 186624S + 3079296C(3) + 4039632 



303264^^^ + 3639168^^ + 325728^' 

V3 



+ 
+ 



176256V''(i) - 1175047r2 



V3 
3494016 



CpTpNf 



138240V'' (Dvr^ - 103680(V''(5))^ + 1537056V''(|) 

- 34992V/"(i) - 1679616s2(f ) + 3359232s2(f ) 
+ 2799360s3(|) - 2239488.S3(f ) + 472327r'* 



- 10247047r^ + 139968S 



1399680C(3) + 729648 

3 



11664^^ + 139968^ + 12528^ 



cl 



72^^^ (a, 0) 



RI'/SMOM 



157464 
(4.14) 

(4.15) 



Again the final expression is the same as that for the vector operator. In numerical form we 

have 



7iT^(a,0) 



7i^^(a,0) 



RI'/SMOM 



RI'/SMOM 



= 3.5555556a + [104.7024244 - 6.5770518iV/] 



+ 



4010.9803829 - 624.8817671iV/ + 14.9653337iV/l + 0{a^) 



= - 1.7777778a - [46.3028612 - 2.7468825A(/-] , 
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- 1692.0143513 - 265.3339715iV/ + 6.0846171Af ] + 0{a'^) 
72'^^K0)^^,,,,_, = Oia') (4.16) 



RI'/SMOM 



for SU{3). 



5 Discussion. 

We have provided the complete set of ampUtudes at two loops for the second moment of the 
Wilson operators inserted in a quark 2-point function in both the MS and RI'/SMOM renor- 
malization schemes. This is not as straightforward a task in comparison with earlier work to 
this level, [28\ [29] , as there is mixing with a total derivative operator. One of the aims of the 
original RI'/SMOM scheme was that the convergence of the conversion functions between these 
two schemes would improve compared with the RI' scheme. Indeed for the quark currents that 
appears to be the case. However, for W2 if anything the RI' result seems to be converging 
marginally quicker. Though it is not completely clear if this is really an appropriate compari- 
son. This is because the operator mixing simply does not arise in the RI' case due to the very 
nature of the momentum configuration used for the Green's function. Indeed given the lack of 
multiplicative renormalizability and hence mixing, it is not entirely clear what the status of the 
RI' scheme is for W2. It may be that RI'/SMOM is the only proper scheme to use of the two. 
Moreover, RI'/SMOM should not suffer the infrared issue associated with RI', [27]. However, 
if one was concerned about improving the convergence of the conversion function it might be 
possible to exploit the freedom one has in actually defining the scheme. For the operator W2 
we used the channel 1 and 2 amplitudes as the basis for the renormalization conditions where 
these channels are defined with respect to a choice of basis tensors. This basis is by no means 
unique as one could equally choose another basis and hence use the analogous channels 1 and 2 
there. As a variation one could instead use a different combination of amplitudes which equates 
to projecting by a linear combination of basis tensors. The advantage of this is that one would 
incorporate more information about the operator within the Green's function which is encoded 
in the other amplitudes. Whilst mathematically this is not inequivalent to a different choice of 
basis tensors, it would avoid the tedious reworking of the construction of the projection matrix 
and thereafter the running of the underlying computer algebra programmes. This is one of the 
reasons why we have provided the MS results so that an interested reader has the information 
for whatever scheme definition variation one might conceive. As was noted in [31] the renor- 
malization of the next moment of the set of Wilson operators involves that of W2, because of 
mixing into a total derivative of W2 itself. Therefore, our results for the anomalous dimensions, 
amplitudes and conversion functions for level W2 will provide important checks on the n = 3 
Wilson operator renormalization in RI'/SMOM, 
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Table 1. Coefficients oiCp for one loop MS Wi amplitudes. 
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Table 2. Coefficients of for one loop MS W2 amplitudes continued. 
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Table 3. Coefficients of Cp for one loop MS dW^ amplitudes. 



24 



(21) 

an 


M/2,(21) 
C(l)n 


MS 


M/2,(21) 
^(2)n 


MS 


14/2,(21) 
C(3)n 


MS 


14/2,(21) 
C(4)n 


MS 


H'2,(21) 
^(5)n 


MS 


1 

TT^a 


730/81 



-1937/162 



8/27 






-8 



-308/27 



a 



























7r2 
V^'(l/3) 
V''(l/3)a 


68/243 
-34/81 



-220/243 
110/81 



80/243 
-40/81 



-248/243 
124/81 



32/27 
32/27 




Table 4. Coefficients of CpTpN^ for two loop MS W2 amplitudes. 
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Table 5. Coefficients of CpTpNf for two loop MS W2 amplitudes continued. 
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Table 6. Coefficients of CpTpNf for two loop MS dW2 amplitudes. 



25 



22 


M/2,(22) 




M/2,(22) 




M/2,(22) 




^2,(22) 




M/2,(22) 




all ' 


C(l)n 


MS 


C(2)n 


MS 


C(3)n 


MS 


C(4)n 


MR 


C(5)n 


MR 


1 


-15361/648 


21445/648 


-188/81 


8351/324 


5879/162 




-37/54 


17/54 


-214/81 


67/81 


266/81 




-2/81 


2/27 


-8/243 


32/243 


76/243 


C(3)a 




1/2 


7/2 




-2 




7/3 




6 


Sa 







2/3 


4/9 




8/9 


16/9 


a 


-13/12 


-107/12 


-23/9 


-125/18 


-107/9 


7r2a2 


1/27 


-7/27 


-10/81 


-26/81 


-58/81 




-5/24 




5/3 




5/9 


-55/36 


-49/18 




























































7/9 


-425/243 


7724/729 


-2384/729 


-7774/729 




4/27 


-22/81 


200/243 


-100/243 


-68/81 


C(3) 


73/6 


-77/6 


98/3 


-41/3 




-34 


s 




1/3 




-1 


4/9 


-10/9 




8/9 


S2(vr/6) 




-43 




53 




208 




70 




168 


52(71/6)0; 




9 




-15 




24 




-26 




-72 


S2(7r/2) 




86 




106 


416 




140 




336 


S2(7r/2)a 




-18 




30 




-48 




52 




144 


S3(7r/6) 


215/3 


-265/3 


1040/3 


-350/3 




280 


53(71/6)0 




-15 




25 




-40 


130/3 




120 


S3(7r/2) 


-172/3 


212/3 


-832/3 


280/3 


224 


S3(7r/2)a 




12 




-20 




32 


-104/3 




-96 


V''(l/3) 




7/6 


425/162 


-3862/243 


1192/243 


3887/243 


V''(l/3)a 


37/36 


-17/36 


107/27 


67/54 


-133/27 


■0 (l/3)a'^ 


-1/18 


7/18 


5/27 


13/27 


29/27 


tp (l/6)a 



























V'(l/3)7r2 


8/27 







64/81 


16/81 







(V''(l/3))2 




2/9 







-16/27 


-4/27 







V''"(l/3) 


-5/54 


11/108 


-11/27 


7/54 


17/54 


<'(l/3)a 


1/108 


-1/36 


1/81 


-4/81 


-19/162 


7r3a/\/3 
7rV^/3 


-29/432 


145/1296 


-29/162 


377/1944 


29/54 


1247/3888 


-1537/3888 


377/243 


-1015/1944 


-203/162 


7rln(3)a/\/3 




3/4 


5/4 




-2 


13/6 




6 


7rln(3)/V^ 

7r(ln(3))2a/\/3 


43/12 


-53/12 


52/3 


-35/6 




-14 


1/16 


-5/48 




1/6 


-13/72 




1/2 


7r(ln(3))VV3 


-43/144 


53/144 


-13/9 


35/72 


7/6 



Table 7. Coefficients of CfCa for two loop MS W2 amplitudes. 
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Table 8. Coefficients of CfCa for two loop MS W2 amplitudes continued. 
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2/3 


4/3 


4/3 


4/3 







a 




-10 


-28/3 




-7 


-14/3 











2/9 




4/9 




4/9 




4/9 


2/27 




-15/8 




-2 




3/2 




-1 































































-236/243 


532/243 


-472/243 


-164/27 


-1672/243 




-10/81 


-64/243 


-20/81 


-56/243 


-32/81 


C(3) 




2/3 




-4 


-22/3 


-32/3 


-28/3 


s 




2/3 




8/3 




4/3 







-28/3 


S2(vr/6) 




10 









20 




40 




68 


S2(7r/6)a 




-6 









-12 




-80 




4 


S2(7r/2) 




86 









-40 




-80 




136 


S2(7r/2)Q 




12 









24 




48 




-8 


S3(7r/6) 


-50/3 







-100/3 


-200/3 


-340/3 


S3(7r/6)a 




10 









20 




40 


-20/3 


S3(7r/2) 


40/3 







80/3 


160/3 


272/3 


S3(7r/2)a 




-8 









-16 




-32 


16/3 


^'{1/3) 


118/81 


-266/81 


236/81 


82/9 


836/81 


ij'{l/3)a 




5/9 


28/9 


10/9 




8/9 




1/3 






1/3 




2/3 




2/3 




2/3 




1/9 


ip [l/6)a 



























V^'(l/3)7r2 


8/27 


64/81 


16/27 


32/81 


16/27 


(V'(l/3))2 




2/9 


-16/27 




4/9 


-8/27 




4/9 


<'(l/3) 


1/108 







1/54 


1/27 


2/27 


V'"(l/3)a 


-1/54 


-1/27 


-1/27 


-4/81 


1/54 


7r3a/\/3 
v^V^/3 


29/648 







29/324 


29/162 


-29/972 


-145/1944 







-145/972 


-145/486 


-493/972 


TT ln(3)a/ \/3 




1/2 









1 




2 




1/3 


7rln(3)/V^ 

7r(ln(3))2a/^/3 




5/6 









5/3 


-10/3 


-17/3 


-1/24 







-1/12 




1/6 


1/36 


7r(ln(3))VV3 


5/72 







5/36 


5/18 


17/36 



Table 9. Coefficients of CfCa for two loop MS dW2 amplitudes. 



28 



23 


IV2,(23) 








M/2,(23) 




IV2,(23) 




IV2,(23) 




ah. 


C(l)n 


MS 


C(2)n 


MS 


C(3)n 


MS 




MS 


''(5) n 


MS 


1 


3971/324 


-9805/648 


1238/81 


-1669/162 


-1987/81 


TT^a 


248/243 


-824/243 


1288/729 


-2296/729 


-6152/729 


TT^a 







16/81 


32/243 


64/243 


128/243 


C(3)a 




4/3 




8/3 




8/3 


8/3 




8 


Sa 







4/3 


8/9 


16/9 


32/9 


a 


-107/108 


539/108 


32/81 


-1697/162 


1838/81 




-2/27 




2/9 


-20/81 


-28/81 


-44/81 































-5/12 


-7/12 


-10/9 


-25/18 


-16/9 




-869/243 


1967/243 


-8852/729 


9362/729 


15424/729 




-68/243 


8/243 


-832/729 


280/729 


608/729 


C(3) 


-46/3 




6 


-136/3 




12 




32 


s 


4/9 


-10/9 


8/27 


-32/27 


-40/27 


S2(7r/6) 




104 




-96 


320 




144 




288 


S2(7r/6)a 




-16 




32 




-32 




32 




96 


S2(vr/2) 




208 




192 




640 


288 




576 


S2(7r/2)Q; 




32 




-64 




64 




-64 




192 


S3(7r/6) 


-520/3 




160 


-1600/3 


240 


480 


S3(7r/6)a 


80/3 


-160/3 


-1600/3 


-160/3 




160 


S3(vr/2) 


416/3 




128 


1280/3 




192 




384 


S3(7r/2)Q; 
^'(1/3) 
V''(l/3)7r2Q 


-64/3 


128/3 


-128/3 


128/3 




128 


869/162 


-1967/162 


4426/243 


-4681/243 


-7712/243 



























V''(l/3)a 


-124/81 


412/81 


-644/243 


1148/243 


3076/243 


V''(l/3)7r^a^ 
V''(1/3)q:^ 





























1/9 




1/3 


10/27 


14/27 


22/27 


V^'(l/3)7r2 


-16/27 







-128/81 


-32/81 


-76/243 


/ / // -I /o\\9 

(1/3))^ 


4/9 







32/27 


8/27 







W (1/3)) a 
(V'(l/3))2a2 




















































<'(l/3) 


29/162 


-1/81 


-23/243 


-76/243 


-76/243 


V'"'(l/3)a 







-2/27 


-4/81 


-8/81 


-16/81 




29/162 


-58/243 


58/243 


-58/243 


-58/81 


-377/486 


58/81 


-580/243 


29/27 


-58/27 


7rln(3)a/\/3 


4/3 




8/3 


8/3 




8/3 




-8 


7rln(3)/V^ 


-26/3 




8 


-80/3 




12 




24 


7r(lii(3))2Q/V3 




1/9 




2/9 




2/9 


2/9 


2/3 


7r(ln(3))VV3 


13/18 




2/3 


20/9 




-1 




-2 



Table 10. Coefficients of Cp for two loop MS W2 amplitudes. 



29 



23 


M/2,(23) 








M/2,(23) 




M/2,(23) 




M/2,(23) 




ah. 


^(6)n 


MS 


C(7)n 


MS 


C(8)r, 


MS 


C(9)„ 


MS 


''(lO) n 


MS 


1 


1609/81 


535/162 


-1994/81 


185/54 


-185/54 


TT^a 


320/729 


-1160/729 


-2368/729 












TT^a 


-32/243 


32/243 


64/243 












C(3)a 




8/3 







8/3 












Sa 




8/9 




8/9 


16/9 












a 


-1406/81 


-401/162 


832/81 




5/6 


5/6 




-4/81 


-20/81 


-28/81 









































4/9 


-11/18 


-14/9 




1/6 




1/6 




-7000/729 


-2774/729 


13604/729 


1858/243 


2786/243 


/I 

vr^ 


-560/729 


-640/729 


64/729 


104/243 


-32/243 


C(3) 


-112/3 


-92/3 


40/3 


28/2 




4/3 


S 


40/27 


-4/27 


-80/27 









4/3 


S2(7r/6) 


224 




160 




224 




-80 




-64 


S2(7r/6)a 




-32 









32 












S2(vr/2) 




448 




320 


448 




160 




128 


S2(7r/2)Q; 




64 









-64 












S3(7r/6) 


-1120/3 


-800/3 


1120/3 


-320/3 


320/3 


S3(7r/6)a 


160/3 







-160/3 












S3(vr/2) 


896/3 


640/3 


-896/3 


-320/3 


-256/3 


S3(7r/2)Q; 
^'(1/3) 
V''(l/3)7r2Q 


-128/3 







128/3 












3500/243 


1387/243 


-6802/243 


-929/81 


-1393/81 



























V''(l/3)a 


-160/243 


580/243 


1184/243 












V''(l/3)7r^a^ 
V''(1/3)q:^ 



























2/27 


10/27 


14/27 







2/9 


V^'(l/3)7r2 


-64/81 


-64/81 







-32/81 


-64/81 


/ / // -I /o\\9 

(1/3))^ 


16/27 


16/27 







8/27 


16/27 


W (1/3)) a 
(V'(l/3))2a2 




















































<'(l/3) 


94/243 


104/243 


-8/243 




1/9 


4/27 


V'"'(l/3)a 


4/81 


-4/81 


-8/81 














58/243 







-58/243 












-406/243 


-290/243 


406/243 


145/243 


116/243 


7rln(3)a/\/3 


8/3 









8/3 












7rln(3)/V^ 


-56/3 


-40/3 


56/3 


20/3 


16/3 


7r(lii(3))2Q/V3 




2/9 









2/9 












7r(ln(3))VV3 


14/9 


10/9 


-14/9 




5/9 




4/9 



Table 11. Coefficients of Cp for two loop MS W2 amplitudes continued. 



30 



(23) 






aiV2,(23) 




yi1/2,(23) 








yW2,(23) 




ah 


C(l)n 


MS 


C(3)n 


MS 


C(4)n 


MS 




MS 


C(9)n 


MS 


1 


-23/8 


-28/3 




-7 


-14/3 









-64/27 


-40/27 


-128/27 




-8 


-4/27 


TT^a 


16/81 


32/81 


32/81 


32/81 







C(3)a 


4/3 







8/3 


16/3 







Ea 


4/3 




8/3 


8/3 




8/3 







a 




4 


32/3 


32/81 


16/3 









-8/27 


-16/27 


-16/27 


-16/27 






































-1 




8/3 




-2 




4/3 









122/27 


176/27 


244/27 


104/9 


172/9 




-20/81 


-256/243 


-40/81 


16/243 


8/27 


C(3) 


-28/3 




-32 


-56/3 


-16/3 




8 


S 




2/3 




8/3 




4/3 









4/3 


S2(7r/6) 




8 




96 




16 




-64 




144 


S2(7r/6)a 




16 









32 




64 







S2(vr/2) 




-16 




192 




-64 




128 


288 


S2(7r/2)a 




-32 







-80/3 




128 







S3(7r/6) 


-40/3 




160 


-160/3 


3203 


240 


S3(7r/6)Q 


-8 


0/3 







64/3 


-320/3 







S3(vr/2) 


32/3 




128 


128/3 


-256/3 




192 


S3(7r/2)a 
V''(l/3) 


64/3 







-121/3 


256/3 







-61/9 


-88/9 


-122/9 


-52/3 


-86/3 


V''(l/3)7r2a 





























32/9 


20/9 


64/9 




12 







it/-. / ^ \ o o 

V''(l/3)7r2Q;2 
^'(1/3)q;^ 



























4/9 




8/9 


8/9 




8/9 


2/9 


V''(l/3)vr2 


-16/27 


-32/27 


-128/81 


-64/81 


-32/27 


(V' (1/3))^ 


4/9 


32/27 


8/9 


16/27 


8/9 


(1/3)) a 
(^'(l/3))2a2 




















































<'(l/3) 




1/6 


32/27 




1/3 




2/7 


1/27 


V^'"(l/3)a 


-2/27 


-4/27 


-4/27 


-8/81 









-29/243 







-58/243 


-116/243 







-29/486 


-58/81 


-29/243 


116/243 


29/27 


7rln(3)a/\/3 




4/3 









8/3 


-16/3 







7rln(3)/^/3 
7r(ln(3))2a/V3 




2/3 




-8 




4/3 


16/3 




12 




1/9 









2/9 




4/9 







7r(ln(3))VV3 


1/18 




2/3 




1/9 




4/9 




-1 



Table 12. CoefHcients of C'p for two loop MS dW2 amplitudes. 



31 



1 


W9 (1) 

^(3)n 


Wo (1) 


W9 (1) 

^(5)n 


Wo (1) 
^(6) n 




Wj (1) 

^(8)n 


Wo (1) 
^(9)n 


Wo (1) 
^(10) n 


1 


16/27 


71/27 


100/27 


-28/27 


37/27 


128/27 


-1/3 


1/3 


2 


-16/81 


-32/81 


-64/81 


16/81 


-16/81 


-32/81 








a 


-8/9 


-16/9 


-26/9 


14/9 


-2/9 


-16/9 








7r2 


-16/243 


64/243 


80/243 


-80/243 


8/243 


160/243 





8/27 


^'(1/3) 


8/81 


-32/81 


-40/81 


40/81 


-4/81 


-80/81 








^'(l/3)a 


8/27 


16/27 


40/27 


-8/27 


8/27 


16/27 





-4/9 



Table 13. Coefficients of Cp for one loop RI'/SMOM W2 amplitudes. 





C(l)n 


C(3) n 


'^(4)n 


C(5)n 


dW2,(l} 
^(9) n 


1 


2 


16/3 


4 


8/3 





TT^a 


-8/27 


-16/27 


-16/27 


-16/27 





a 


-1 


-8/3 


-2 


-4/3 





7r2 


4/27 


16/27 


8/27 





8/27 


^'(1/3) 


-2/9 


-8/9 


-4/9 





-4/9 




4/9 


8/9 


8/9 


8/9 






Table 14. Coefficients of Cp for one loop RI'/SMOM dW2 amplitudes. 



(21) 

ail 


M/2,(21) 
^(3)n 


M/2,(21) 


H'2,(21) 
^(5)n 


M/2,(21) 
^(6)n 


1 


8/27 


-8 


-308/27 


76/27 




320/729 


640/729 


1280/729 


-320/729 


a 


160/81 


320/81 


520/81 


-280/81 




80/243 


-248/243 


-64/81 


64/81 


t//(l/3) 


-40/81 


124/81 


32/27 


-32/27 


(/■'(l/3)a 


-1(30/243 


-320/243 


-(340/243 


1(30/243 



Table 15. Coefficients of CpTrNf for two loop RI'/SMOM W2 amplitudes. 



32 



21) 


^(7) n 




WoA21) 
yy z^y-^^i 

^(9)n 


Wo, (21) 
^(10) n 


1 


-44/9 


-472/27 


32/27 


-32/27 


TT^a 


320/729 


640/729 








a 


40/81 


320/81 








7r2 


-56/243 


-688/243 


8/81 


-232/243 


V^'(l/3) 


28/81 


344/81 


-4/27 


116/81 


V''(l/3)a 


-160/243 


-320/243 









Table 16. Coefficients of CpTpN} for two loop RI'/SMOM W2 amplitudes continued. 



21 

On 


C(l)n 


C(3)n 


C>IV2,(21) 

C(4)n 




C(9)n 


1 


-58/9 


-464/27 


-116/9 


-232/27 







160/243 


320/243 


320/243 


320/243 





Q 


20/9 


160/27 


40/9 


80/27 







-152/243 


-608/243 


-304/243 





-208/243 


^^'(1/3) 


76/81 


304/81 


152/81 





104/81 


V''(l/3)a 


-80/81 


-160/81 


-160/81 


-160/81 






Table 17. Coefficients of CfTeNj for two loop RI'/SMOM dW2 amplitudes. 



33 



22) 


(22) 

^(3)n 


W-> (22) 


Wo (22) 


Wo (22) 


1 


-188/81 


8351/324 


5879/162 


-1637/162 




-2314/729 


-173/729 


842/729 


-1574/729 


4 

TT a 


-8/243 


32/243 


76/243 


-52/243 


C(3)« 


-2 


7/3 


6 


-6 


Sa 


4/9 


8/9 


16/9 


-4/9 


a 


-401/81 


-1901/162 


-3187/162 


1849/162 




-2/9 


-14/27 


-10/9 


10/27 




-1 


-29/12 


-25/6 


5/2 




-4/81 


-8/81 


-16/81 


4/81 


a3 


-2/9 


-4/9 


-13/18 


7/18 




7724/729 


-2384/729 


-7774/729 


3346/729 




200/243 


-100/243 


-68/81 


148/243 


C(3) 


98/3 


-41/3 


-34 


70/3 


s 


4/9 


-10/9 


-8/9 


8/9 


S2(7r/6) 


-208 


70 


168 


-128 


S2(7r/6)a 


24 


-26 


-72 


48 


S2(7r/2) 


416 


-140 


-336 


256 


S2(vr/2)a 


-48 


52 


144 


-96 


S3(7r/6) 


1040/3 


-350/3 


-280 


640/3 


S3(7r/6)a 


-40 


130/3 


120 


-80 


S3(7r/2) 


-832/3 


280/3 


224 


-512/3 


S3(7r/2)a 


32 


-104/3 


-96 


64 


V''(l/3) 


-3862/243 


1192/243 


3887/243 


-1673/243 


V''(l/3)a 
V''(l/3)a2 


1157/243 


173/486 


-421/243 


787/243 


1/3 


7/9 


5/3 


-5/9 


-0 (1/3)q'^ 


2/27 


4/27 


8/27 


-2/27 


^'(1/3)^2 


64/81 


16/81 


17/54 


32/81 


(V^'(l/3))2 


-16/27 


-4/27 


-19/162 


-8/27 


V''"(l/3) 


-11/27 


7/54 


17/54 


-5/18 


V''"(l/3)a 


1/81 


-4/81 


-19/162 


13/162 




-29/162 


377/1944 


29/54 


-29/81 


377/243 


1015/1944 


-203/162 


232/243 


7rln(3)a/\/3 


-2 


13/6 


6 


-4 


7rln(3)/\/3 

7r(ln(3))2a/\/3 


52/3 


-35/6 


-14 


32/3 


1/6 


-13/72 


-1/2 


1/3 


7r(ln(3))VV3 


-13/9 


35/72 


7/6 


-8/9 



Table 18. Coefficients of CfCa for two loop RI'/SMOM W2 amplitudes. 



34 



(22 


W'? (22) 


Wo f22) 


Wo (22) 


Wo i22) 










''(10) n 


1 


4375/324 


4430/81 


-379/108 


379/108 


TT^a 


-1531/729 


-362/729 


7/27 


-1/27 


4 

TT a 


-8/243 


32/243 





-4/81 


C(3)a 


-7/3 


2 


1/3 


-1 


Sa 


4/9 


8/9 








a 


-53/81 


-937/81 


-1/6 


1/6 




-2/9 


-14/27 





2/27 




-1/12 


-7/3 


-1/12 


1/12 




-4/81 


-8/81 










-1/18 


-4/9 










968/729 


-6128/729 


-284/81 


-820/243 




40/243 


-88/81 


-44/243 


-52/243 


C(3) 


19/3 


-110/3 


-5 


-13/3 


S 


-2/9 


-28/9 


-2/9 


-10/9 


S2(7r/6) 


-50 


208 


30 


38 


S2(7r/6)a 


14 


-24 


-2 


6 


S2(7r/2) 


100 


-416 


-60 


-76 


S2(vr/2)a 


-28 


48 


4 


-12 


S3(7r/6) 


250/3 


-1040/3 


-50 


-190/3 


S3(7r/6)a 


-70/3 


40 


-50 


-10 


S3(7r/2) 


-200/3 


832/3 


10/3 


152/3 


S3(7r/2)a 


56/3 


-32 


40 


8 


^'(1/3) 


-484/243 


3064/243 


142/27 


410/82 


V''(l/3)a 


1531/486 


181/243 


-7/18 


1/18 


V''(l/3)a2 


1/3 


7/9 





-1/9 


V' (l/3)a'^ 


2/27 


4/27 








V''(l/3)7r2 


32/81 





16/81 


32/81 


(V''(l/3))2 


-8/27 





-4/27 


-8/27 


<'(l/3) 


-1/9 


11/27 


7/162 


5/162 


V'"'(l/3)a 


1/81 


-4/81 





1/54 




-203/1944 


29/162 


29/1944 


-29/648 


725/1944 


-377/243 


-145/648 


-551/648 


TT ln(3)a/ \/3 


-7/6 


2 


1/6 


-1/2 


7rln(3)/V^ 
7r(ln(3))2Q/V3 


25/6 


-52/3 


-5/2 


-19/6 


7/72 


-1/6 


-1/72 


1/24 


7r(ln(3))VV3 


-25/72 


13/9 


5/24 


19/72 



Table 19. Coefficients of CfCa for two loop RI'/SMOM W2 amplitudes continued. 
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(22 


(22) 

^(l)n 


dW-? (22) 
^(3)n 


dW-> (22) 


dW-> (22) 


dW-> (22) 


1 


707/36 


1414/27 


707/18 


707/27 





TT^a 


-284/243 


-892/243 


-568/243 


-244/243 


2/9 


4 

TT a 


4/81 


8/81 


8/81 


8/81 


-4/81 


C(3)a 














-2/3 


Sa 


2/3 


4/3 


4/3 


4/3 





a 


-223/36 


-446/27 


-223/18 


-223/27 







-10/27 


-20/27 


-20/27 


-20/27 


2/27 




-5/4 


-10/3 


-5/2 


-5/3 







-2/27 


-4/27 


-4/27 


-4/27 







-1/4 


-2/3 


-1/2 


-1/3 







-236/243 


532/243 


-472/243 


-164/27 


-1672/243 




-10/81 


-64/243 


-20/81 


-56/243 


-32/81 


C(3) 


-11/3 


-4 


-22/3 


-32/3 


-28/3 


s 


-2/3 


-8/3 


-4/3 





-4/3 


S2(vr/6) 


10 





20 


40 


68 


S2(7r/6)Q; 


-6 





-12 


-24 


4 


S2(7r/2) 


-20 





-40 


-80 


-136 


S2(7r/2)a 


12 





24 


48 


-8 


S3(7r/6) 


-50/3 





-100/3 


-200/3 


-340/3 


S3(7r/6)a 


10 





20 


40 


-20/3 


S3(7r/2) 


40/3 





80/3 


160/3 


272/3 


S3(7r/2)Q! 


-8 





-16 


-32 


16/3 


V''(l/3) 


118/81 


-266/81 


236/81 


82/9 


836/81 


V''(l/3)a 


142/81 


446/81 


284/81 


122/81 


-1/3 


■0'(l/3)a2 


5/9 


10/9 


10/9 


10/9 


-1/9 


(l/3)a'^ 


1/9 


2/9 


2/9 


2/9 





V''(l/3)vr2 


-8/27 


64/81 


16/27 


32/81 


16/27 


(V''(l/3))2 


-2/9 


-16/27 


-4/9 


-8/27 


-4/9 


<'(l/3) 


1/108 





1/54 


1/27 


2/27 


V''"(l/3)a 


-1/54 


-1/27 


-1/27 


-1/27 


1/54 




29/648 





29/324 


29/162 


-29/972 


-145/1944 





-145/972 


-145/486 


-493/972 


7rln(3)a/\/3 


1/2 





1 


2 


-1/3 


7rln(3)/v/3 

7r(ln(3))2a/V3 


-5/6 





-5/3 


-10/3 


-17/3 


-1/24 





-1/12 


-1/6 


1/36 


7r(ln(3))VV3 


5/72 





5/36 


5/18 


17/36 



Table 20. Coefficients of CfCa for two loop RI'/SMOM dW2 amplitudes. 
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23 


W2,(23) 








' 


^(3)n 






^(6)n 


1 


718/81 


-443/162 


-683/81 


737/81 




1024/729 


-4336/729 


-10376/729 


2960/6561 


4 

71 a 


928/6561 


704/6561 


1404/6561 


608/6561 


C(3)a 


-8/3 


8/3 


8 


-8/3 


Sa 


8/9 


16/9 


32/9 


-8/9 


q; 


-188/81 


-125/162 


238/81 


-454/81 




116/243 


220/243 


380/243 


-140/243 




128/2187 


256/2187 


512/2187 


-128/2187 




26/27 


89/54 


68/27 


-32/27 




-30452/2187 


31286/2187 


52576/2187 


-26440/2187 




-24064/19683 


8776/19683 


18656/19683 


-17360/19683 


C(3) 


-136/3 


12 


32 


-112/3 


s 


8/27 


-32/27 


-40/27 


40/27 


S2(7r/6) 


320 


-144 


-288 


224 


S2(vr/6)a 


-32 


32 


96 


-32 


S2(vr/2) 


-640 


288 


576 


-448 


S2(7r/2)a 


64 


-64 


-192 


64 


S3(7r/6) 


-1600/3 


240 


480 


-1120/3 


S3(7r/6)a 


-1600/3 


-160/3 


-160 


160/3 


S3(7r/2) 


1280/3 


-192 


-384 


896/3 


S3(7r/2)a 


-128/3 


128/3 


128 


-128/3 


^'(1/3) 

'(/''(l/3)7r2a 


15226/729 


-15643/729 


-26288/729 


13220/729 


-64/2187 


1024/2187 


2048/2187 


-1472/2187 


V''(l/3)a 


-512/243 


2186/243 


5188/143 


-1480/243 


V''(l/3)7r2a2 


-128/729 


-256/729 


-512/729 


128/729 


V''(l/3)a2 


-58/81 


-110/81 


-190/81 


70/81 


V''(l/3)^2 


-8768/6561 


-3808/6561 


-2240/6561 


-2944/6561 


(^'(l/3))2 


2192/2187 


952/2187 


560/2187 


736/2187 


/ ///I /o^^2 
(V' (1/3)) a 

(V'(l/3))2a2 


lb/1 IS) 


— 25d/7z9 


—512/729 


o/^o /'von 

368/729 


32/243 


64/243 


128/243 


-32/243 


V''"(l/3) 


152/243 


-23/243 


-76/243 


94/243 


<'(l/3)a 


-4/81 


-8/81 


-16/81 


4/81 




58/243 


-58/243 


-58/81 


58/243 


-580/243 


29/27 


-58/27 


-406/243 


7rln(3)a/\/3 


8/3 


-8/3 


-8 


8/3 


7rln(3)/V^ 
7r(ln(3))2a/\/3 


-80/3 


12 


24 


-56/3 


-2/9 


2/9 


2/3 


-2/9 


7r(ln(3))V^/3 


20/9 


-1 


-2 


14/9 



Table 21. Coefficients of for two loop RI'/SMOM W2 amplitudes. 
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23 


W2,(23) 




W2,(23) 


1^2,(23) 


ah. ' 




^(8)n 


^(9)n 


^(10) n 


1 


605/162 


-610/81 


83/54 


-83/54 




-1496/729 


-5272/729 





-16/27 


4 

IT a 


736/6561 


-64/6561 





-64/729 


C(3)a 





8/3 










8/9 


16/9 








a 


-523/162 


-244/81 


-7/18 


7/18 




68/243 


220/243 





-4/27 




128/2187 


256/2187 








a2 


19/54 


46/27 


1/6 


-1/6 




-8930/2187 


49028/2187 


566/81 


1030/81 




-17632/19683 


5056/19683 


872/2187 


-128/2187 


C(3) 


-92/3 


40/3 


28/3 


-4/3 




-4/27 


-80/27 





-4/3 


S2(7r/6) 


160 


-224 


-80 


128 


S2(7r/6)a 





32 








S2(vr/2) 


-320 


448 


160 


320/3 


52(71/2)0 





-64 








S3(7r/6) 


-800/3 


1120/3 


400/3 


-256/3 


S3(7r/6)a 





-160/3 








S3(7r/2) 


896/3 


-896/3 


-320/3 


-515/27 


53(71/2)0; 





128/3 








V''(l/3) 


4465/729 


-24514/729 


-283/27 


-515/27 


V''(l/3)7r2o 


128/2187 


1792/2187 





64/243 


^'(l/3)a 


748/243 


2636/243 





8/9 


V''(l/3)7r2o2 


-128/729 


-256/729 








^'(1/3)0^ 


-34/81 


-110/81 





2/9 


V''(l/3)7r2 


-4821/6561 


-3328/66561 


-224/729 


-736/729 


(V''(l/3))' 


1208/2187 


832/2187 


56/243 


184/243 


(V' (1/3))^q; 
(^'(l/3))2o2 


—32/729 


A AO /'7or» 

—448/729 





— 16/81 


32/243 


64/243 










104/243 


-8/243 


-1/9 


4/27 


V'"'(l/3)o 


-4/81 


-8/81 








7r3o/\/3 
7rVx/3 





-58/243 








-290/243 


406/243 


145/243 


116/243 


7r ln(3)a/ \/3 





-8/3 








7rln(3)/^/3 
7r(ln(3))2a/^/3 


-40/3 


56/3 


20/3 


16/3 





2/9 








7r(ln(3))VV3 


10/9 


-14/9 


-5/9 


-4/9 



Table 22. Coefficients of Cj. for two loop RI'/SMOM W2 amplitudes continued. 
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23 








dWo (231 


dWo (231 




^(l)n 


^(3)n 


^(4)n 






1 


-7/2 


-28/3 


-7 


-14/3 







-68/27 


-56/27 


-136/27 


-8 


-8/27 


4 

TT a 


16/81 


32/81 


32/81 


32/81 





C(3)a 


4/3 














Eq 


4/3 


8/3 


8/3 


8/3 





a 


2 


16/3 


4 


8/3 



















-4/27 






































122/27 


176/27 


244/27 


104/9 


172/9 




-20/81 


-256/243 


-40/81 


16/243 


8/27 


C(3) 


-28/3 


-32 


-56/3 


-16/3 


8 


s 


-2/3 


-8/3 


-4/3 





-4/3 


S2(7r/6) 


8 


96 


16 


-64 


-144 


S2(7r/6)Q! 


16 





32 


64 





S2(vr/2) 


-16 


-192 


-32 


128 


288 


S2(7r/2)a 


-32 





-64 


-128 





S3(7r/6) 


-40/3 


-160 


-80/3 


320/3 


240 


S3(7r/6)Q! 


-80/3 





-160/3 


-320/3 





S3(vr/2) 


32/3 


128 


64/3 


-256/3 


192 


S3(7r/2)Q! 


64/3 





128/3 


256/3 





V''(l/3) 

V''(l/3)7r2a 


-61/9 


-88/9 


-122/9 


-52/3 


-86/3 

















V''(l/3)a 


34/9 


28/9 


68/9 


12 


4/9 


V''(l/3)7r2a2 

















^'(l/3)a2 














-32/27 


V''(l/3)vr2 


-16/27 


-128/81 


-32/27 


-64/81 


-32/27 


(V''(l/3))2 


4/9 


32/27 


8/9 


16/27 


8/9 


(V' (1/3)) a 
(^'(l/3))2a2 


— 2/27 


—A/n 


U 











32/27 


8/9 








vy"(i/3) 


4/9 


16/27 


1/3 


2/27 


1/27 


V'"'(l/3)a 


-2/27 


-4/27 


-4/27 


-4/27 





TT^a/V^ 
7rVx/3 


-29/243 





-58/243 


-116/243 





-29/486 


-58/81 


-29/243 


116/243 


29/27 


7rln(3)a/\/3 


-4/3 





-8/3 


-16/3 





7rln(3)/^/3 
7r(ln(3))2a/^ 


-2/3 


-8 


-4/3 


16/3 


12 


1/9 





2/9 


4/9 





7r(ln(3))VV3 


1/18 


2/3 


1/9 


-4/9 


-1 



Table 23. Coefficients of C|, for two loop RI'/SMOM 51^2 amplitudes. 
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